Terms in elliptic divisibility sequences divisible by their 

indices 
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Abstract. Let D = (D n ) n >i be an elliptic divisibility sequence. 
We study the set <S(D) of indices n satisfying n | D n . In par- 
ticular, given an index n £ 5(D), we explain how to construct 
elements nd £ 5(D), where d is either a prime divisor of D n , or d 
is the product of the primes in an aliquot cycle for D. We also give 
bounds for the exceptional indices that are not constructed in this 
way. 



Introduction 

In this note we investigate the terms in elliptic divisibility sequences 
that are divisible by their indices. The analogous problem has been 
studied for a number of other types of sequences. For example, the 
Fibonacci sequence (F n ) n >i satisfies 

n | F n n G {1,5, 12, 24, 25, 36, 48, 60, 72, 96, . . .}. 

See [1, 9, 10, 19, 20, 21, 24] for results on index divisibility in the 
Fibonacci sequence and in more general Lucas sequences. To cite an- 
other example, values of n that divide a n — a are called pseudoprimes 
to the base a. They have been studied for their intrinsic interest and 
for applications to cryptography [2, 11, 12, 14, 23]. 

In general, for any integer sequence A = (A n ) n >i we define the index 
divisibility set of A to be 

«S(A) = {n>l:n\ A n ). 

Our goal is to build 5(A) multiplicatively via a directed graph that 
connects each element n G 5(A) to its (minimal) multiples in 5(A). 
Thus we define a directed graph by taking the set 5(A) to be the set 
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of vertices and by drawing an arrow from n to to if the following two 
conditions are true: 

(1) n | to. 

(2) If k e S(k) satisfies n \ k \ m, then k = n or k = to. 

In other words, if we partially order S(A) by divisibility, then we draw 
an arrow from n to to if n is strictly smaller than to and if there are 
no elements of S(A) that are strictly between n and to. 

We denote the set of arrows by Arr(A), and we assign weight m/n 
to the arrow (n — > to). (Smyth [19, Section 8] defines a similar struc- 
ture, but he allows only arrows of prime weight, so his graphs may be 
disconnected.) 

Definition. Let E/Q be an elliptic curve given by a Weierstrass equa- 
tion and let P G E(Q) be a nontorsion point. The elliptic divisibility 
sequence (EDS) associated to the pair (E, P) is the sequence of positive 
integers D = (D„) n >i obtained by writing 



as a fraction in lowest terms. The EDS is minimal if E is given by 
a minimal Weierstrass equation. An EDS is normalized if D 1 = 1. 
An arbitrary EDS (D n ) n >i can be normalized by a change of vari- 
ables in the defining Weierstrass equation, in which case the new EDS 
is (D n / Di) n >i. Note, however, that the normalized sequence may not 
be minimal. 

We remark that there is an alternative definition of EDS via a non- 
linear recurrence that gives almost the same set of sequences; see Re- 
mark 28 for further details. We also note that, as its name suggests, 
an EDS is a divisibility sequence, i.e., 



The arithmetic properties of EDS have been extensively studied as 
examples of nontrivial nonlinear recursions that possess enough addi- 
tional structure to make them amenable to Diophantine analysis. See 
for example Ward's original papers [25, 26], subsequent work includ- 
ing [5, 7, 8, 18], and applications of EDS to Hilbert's 10th problem and 
to cryptography [3, 6, 13, 22]. 

Although EDS are defined via a non-linear process, their underlying 
structure comes from the associated elliptic curve. They are thus a 
natural generalization of linear recursions such as the Fibonacci and 
Lucas sequences, which are associated to the multiplicative group. 
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Example 1. Let D be the EDS 

D = (1, 1, 1, 1, 2, 1, 3, 5, 7, 4, 23, 29, 59, 129, 

314, 65, 1529, 3689, 8209, 16264, 83313, . . . ) 
associated to the elliptic curve and point 

E : y 2 + y = x 3 - x, P = (0,0). 

Then 

«S(D) = {1, 40, 53, 63, 80, 127, 160, 189, 200, 320, 400, 441, 443, . . . }. 

We remark that the sequence D grows very rapidly. Thus the first two 
nontrivial elements of S(D) in this example come from 

£>4o = 40 • 13526278251270010, 

D 53 = 53 • 299741133691576877400370757471. 

The reader may have noticed that <S(D) contains the primes 53, 127, 
and 443, which are the first three anomalous primes for E, i.e., primes 
satisfying #E(F P ) = p. This is not a coincidence. 

Smyth has given an explicit description of index divisibility for Lu- 
cas sequences. For comparison with our results, we state one of his 
theorems, reformulated using the terminology of directed graphs. 

Theorem 2. (Smyth [19, Theorem 1]) Let a,b E Z, and let L = 

(L n ) n >i be the associated Lucas sequence of the first kind, i.e., defined 
by the recursion 

L n +2 — aL n+ i — bL n , L = 0, L\ — 1. 

Let A = a 2 — 4b. Then the arrows originating at a vertex n e «S(L) are 

{n np : p is prime and p \ L n A} U B a ,b, 

where 

{{1 — y 6} if a = 3 (mod 6) and b = ±1 (mod 6), 
{1 ->■ 12} if a = ±1 (mod 6) and b = -1 (mod 6), 
otherwise. 

Smyth's theorem says in particular that with at most one exception, 
every arrow for a Lucas sequence has prime weight. This is not true 
for EDS and is due to the fact that the number of points #E(W q ) 
on an elliptic curve over a finite field varies irregularly compared to 
the number of points in the multiplicative group F* of a finite field. 
This leads to EDS arrows of the form n —y nd, where d is a so-called 
aliquot number for the EDS, as in the following definition. The aliquot 
phenomenon has no analogue in the case of Lucas sequences. 



4 



JOSEPH H. SILVERMAN AND KATHERINE E. STANCE 



Definition. A list (pi,---,pe) of distinct primes of good reduction 
for E is an aliquot cycle for D if 

Pi+i = min{r > 1 : Pi \ D r } for all 1 < i < £, 

where we set pe + ± = p\ to complete the cycle. The associated aliquot 
number is the product p\ • • - pi. 

The index divisibility graph of an EDS is considerably more compli- 
cated than that of a Lucas sequence. We state here a simplified version 
of Theorem 15, which is the main result of this paper. We remark that 
an analogue of our main result for EDS associated to singular elliptic 
curves would give a version of Smyth's theorem; see Remark 24 for 
details. 

Theorem 3. Let D be a minimal regular EDS associated to the elliptic 
curve E/Q and point P G E(Q). (See Section 2 for the definition of 
regularity. In particular, every EDS has a regular subsequence.) 

(a) If n G S(D) and p is prime and p \ D n , then (n — > np) G Arr(D). 

(b) If n G 5(D) and d is an aliquot number for D and gcd(n, d) = 1, 
then (n — > nd) G Arr(D). 

(c) If p > 7 is a prime of good reduction for E and if (n — > np) G 
Arr(D), then either p | D n or p is an aliquot number for D. 

(d) If gcd(n, d) = 1 and if (n — > nd) G Arr(D) and if d = P1P2 • • -pe is 
a product of t > 2 distinct primes of good reduction for E satisfying 
mm pi > (2 -1 / 2 ^ — l)~ 2 , then d is an aliquot number for D. 

We briefly describe the contents of this note. In Section 1 we give 
some basic properties of elliptic divisibility sequences. In particular, 
Lemma 5 states fairly delicate divisibility estimates whose origins lie in 
the formal group of E. The brief Section 2 gives the definition of aliquot 
cycles and aliquot numbers for EDS. Section 3 contains the statement 
and proof of Theorem 15, which is the main result of this paper. The- 
orem 15, which is an expanded version of Theorem 3, explains how to 
construct the arrows that are used to build 5(D). This is followed in 
Section 4 with a number of remarks and examples related to our main 
theorem. Section 5 defines aliquot cycles on an elliptic curve (see [17]) 
and explains how they are related to aliquot cycles for an EDS on that 
curve. Finally, in Section 6, we make some miscellaneous remarks on 
general index divisibility sets and on an alternative definition of EDS. 

1. Preliminaries on elliptic divisibility sequences 

Let D be a minimal EDS associated to an elliptic curve E/Q and 
point P G E(Q). We let Disc(£ l ) denote the minimal discriminant 
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of E. For all primes p we have 

p | D n •<=>- [n]P = (mod p). 

Definition. We write r n = r n (D) for the rank of apparition of n in D, 
which is defined by 

r n = min{r > 1 : n \ D r }. 

Let <£/ SpecZ denote the Neron model of E. Then an equivalent defi- 
nition of r n is that it is the smallest value of r > 1 such that 

[r]P = O (mod n), 

where the congruence takes place in (£(Z/nZ). 

The following three lemmas contain virtually all of the information 
about EDS that we will use in our analysis of EDS index divisibility. 

Lemma 4. Let D be a minimal EDS associated to an elliptic curve E/Q 
and point P G -E'(Q). Then 

n | D m r n | m <^=^> [m]P = O (mod n). 

Proof. Immediate from the definitions. □ 

The next lemma describes the growth of p-divisibility for EDS. A 
direct corollary is that an EDS is a divisibility sequence. 

Lemma 5. Let D = (-D„) n >i be a minimal EDS, let n > 1, and let p 
be a prime satisfying p \ D n . 

(a) For all m > 1 we have 

ord p (D mn ) > ord p (mD n ). 

(b) The inequality in (a) is strict, 

ord p (D mn ) > ord p (mD n ), 

if and only if 

, W7^\ i , f E has ordinary or multi-\ 
p = 2 and 2 m and ordol-LL) = 1 and , n ■ 

1 v ' y phcative reduction at 2 J 

(For the definition of ordinary reduction, see [16, §V.3]. In partic- 
ular, E has ordinary reduction at 2 if and only if 2 | #E(W2).) 

Proof. The assumption that p \ D n is equivalent to the assertion that 
[n]P is in £ , 1 (Q P ), the kernel of reduction modulo p. We use the stan- 
dard isomorphism between Ei(Q p ) and the formal group E(pZ p ) asso- 
ciated to E given by 

: E 1 (Q p ) — > E(pZ p ), (x,y) — )■ -x/y. 
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Note that this isomorphism is valid even if E has bad reduction at p, 
in which case E is the formal additive or multiplicative group. (See 
[16, Chapter IV] for basic properties of formal groups.) 

Writing x([n]P) = (A n / D 2 , B n / D 3 ), our assumption that p | D n 
implies that p \ A n B n , so 

ord p 0([n]P) = ord p (- A n D n /B n ) = ord p (D n ). (1) 

Standard properties of formal groups [16, IV. 2. 3(a), IV. 4. 4] say that 
the multiplication-by-p map has the form 

iph(z)= P f(z)+g(*), (2) 

where f,g e 7, p [[z\] are power series with no constant term, and / has 
the form f(z) = z + 0(z 2 ). It follows that for ord p (,2) > 1, we have 

ord p ([p]g(z)) > ord p (p4 (3) 

We write m = p k s with p\ s. Repeated application of (3) gives 

ord p (\p%(z))>ord p (p k z). (4) 

Further, we have [s]e(z) = sz + 0(z 2 ), so 

ord p ([sp k ]z(z)) =ord p ([p k ]zz). (5) 

Combining (4) and (5) gives 

ord p ([m]^(z)) > ord p (p k z), 

with equality if k — 0. Substituting z = 0([n]P) and using (1) 
gives (a), and it also gives (b) if p \ m. 

To prove (b) in general, we assume that p \ m. Analyzing (2) more 
closely, we see that 

ord p([pk(^)) = oid p (pz) = oTd p {z) + 1 (6) 

unless ordp(pz) =pord p (z). (Note that ord p (,2) > 1.) Since 

ord p (pz) = pord p (z) 1 = (p - 1) oid p (z), 

we see that (6) holds except possibly in the case p = 2 and ord p (z) = 1. 

Suppose now that p = 2 and ord p (,2) = 1. The formal group law for 
an elliptic curve starts [16, §IV.l] 

[2]g(z) =2z - a x z 2 - 2a 2 z 3 + (a 3 + a 1 a 2 )z 4 H , 

where a\, . . . ,a& are Weierstrass coefficients. Hence under the assump- 
tion that ord 2 (z) > 1, we see that (6) fails if and only if 

ord 2 (2^ - aiz 2 + 0{2z 3 )) > 3 <^ ord 2 (l - a lZ /2) > 1 

<^=^ ord 2 (ai) = 0, i.e., a\ e 1T 2 . 
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(The last implication follows because z = 2 (mod 4), so 1 — a\z/2 = 
1 — ai (mod 2).) If E has good reduction modulo 2, then j(E) = 
aP/Disc(£) (mod 2), so [16, Exer. 5.7] gives 

ord 2 (cii) = j(.E') ^ (mod 2) £ is ordinary mod 2. 

On the other hand, if E has bad reduction at 2, then an easy compu- 
tation shows that a± = 1 (mod 2) for multiplicative reduction and 
ai = (mod 2) for additive reduction. This completes the proof 
that (6) fails if and only if p — 2 and p | m and ord p (-D n ) = 1 and 
E has either ordinary or multiplicative reduction. We call this the 
exceptional case. 

Repeated application of (6) shows that if we are not in the excep- 
tional case, then 

ord p ([p k ]z(z)) = ordp(z) + k. 

In the exceptional case, the first multiplication by [p] gives a strict 
inequality, after which we are out of the exceptional case and can ap- 
ply (6), so we find that 

OYd p ([p k ]^{z)) = ordp([p]^(z)) + k - 1 > ord p (,2) + k. 

Now using (5) and the fact that m = p k s with p \ s, we get 

ordp([m\fi(z)) > ord p (mz) 

in the exceptional case and 

ordp([m\z(z)) = oid p (mz) 

otherwise. Substituting z = 0([n]P) and using (1) proves (b). □ 

The third lemma gives bounds for r p . 

Lemma 6. Let D be a minimal EDS associated to an elliptic curve E/Q 
and point P e -E'(Q) and let p be a prime. Then 

r n | #(£(Z/nZ). 

In particular, if p is a prime with P e E QS (¥ P ), then 

r P < (v^+1) 2 . 

If P G E ns (¥ p ) and E has bad reduction at p, then r p divides p — 1, 
p + 1, or p depending respectively on whether the reduction is split 
multiplicative, non-split multiplicative, or additive. 

Proof. The first statement is immediate, since r n is the order of the 
point P in the group (£(Z/nZ). The estimates for r p follow from the 
Hasse-Weil bound #E(W P ) < (y/p + l) 2 when E has good reduction, 
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and the explicit description of E ns (¥ p ) for the three types of bad re- 
duction. □ 

Example 7. The minimal EDS associated to 

E : y 2 + xy = x 3 - 2x + 1 and P=(1,0) 
is the sequence 

1,1,1,2,1,3,7,8,25,37,.... 
Thus £> 4 = 2 and D s = 8, so 

3 = ord 2 (D 2 ») > ord 2 (2£> 2 2) = 2. 

The strictness of the inequality in Lemma 5(a) corresponds to the 
exceptional case p = 2, to = 2, and n — 4, where we note that 
ord 2 (-D 4 ) = 1 and #P(F 2 ) = 4, so in particular E has ordinary re- 
duction at 2. 

Remark 8. More generally, for any integer N > 2 there exists a mini- 
mal EDS such that 

ord 2 (D 2 ) = ord 2 (Dx) + N. 

Here is one construction. Choose an elliptic curve of positive rank 
having a rational 2-torsion point T in the formal group £ , (2Z 2 ). Taking 
a multiple of a point of infinite order, we can find a rational nontorsion 
point Q in E(2 N Z 2 ). Then the EDS associated to P = Q + T will have 
ord 2 (.Di) = 1 and ord 2 (D 2 ) — N + 1. The reason that this only works 
for the prime p = 2 is because for p > 3, the formal group E(7L P ) is 
torsion free; in fact, it is isomorphic to the additive group Z+. 

Proposition 9. Let D be a minimal EDS. 

(a) D is a divisibility sequence. 

(b) The set S(D) is closed under multiplication. 

Proof, (a) We need to prove that D m \ D rnn . It suffices to prove that 
ord p (D mn ) > ord p (-D n ) for all primes p, but this is immediate from 
Lemma 5(a). 

(b) Suppose that to, n G S(D) and let p \ n. Then p \ D n , so 
Lemma 5(a) and the assumption that n \ D n give 

ordp(-D m „) > ord p (TO.D n ) > ord p (TOn). 

Reversing the roles of to and n for p \ to again gives ord p (-D mn ) > 
ordp(TOn). Hence mn \ D mn , so mn e5(D). □ 

Remark 10. If p > 3 and p \ Di, then Lemma 5(b) with n = 1 says 
that ord p (-D m ) = oid p {mDi) for all to. 
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Remark 11. Although we will not need this fact, we mention that ellip- 
tic divisibility sequences grow extremely rapidly. Thus if D is associated 
to (E,P), then 

hm — = h E (P), 

n— >oo fl 

where h E (P) > is the canonical height of P [16, VIII §9]. 

2. Aliquot Cycles and Aliquot Numbers for EDS 

In this section we define aliquot cycles and aliquot numbers associ- 
ated to an EDS. 

Definition. Let D be an EDS associated to the curve E(Q) and point 
P G E(Q). We recall that r n (D) denotes the rank of apparition of n 
in the sequence D; see Section 1. An aliquot cycle {of length I) for D 
is a sequence (pi,p 2 , • • • ,Pe) of distinct primes of good reduction for E 
such that 

r pi (D)=p 2 , r p2 (D)=p 3 , r^^D) = p e , r pi (D) = p 1 . 

An amicable pair is an aliquot cycle of length two. 

If we drop the requirement that E have good reduction, then we call 
(Pi>P2, • • • ,Pe) a generalized aliquot cycle. 

In our study of index divisibility for EDS, the products of the primes 
appearing in each aliquot cycle play a key role, so we give them a name. 

Definition. Let D be a minimal EDS. We define the set of aliquot 
numbers of D to be 

A(D) = {p 1 ■ ■ -pe : (pi, . . . ,pt) is an aliquot cycle for D}. 

We also define the larger set 

^4 gcn (D) = {pi ■ ■ - pe : (pi, . . . ,pi) is a generalized aliquot cycle for D}. 

Remark 12. We observe that an aliquot cycle of length one consists of 
a single prime p satisfying r p (D) = p. If p > 7, Hasse's estimate for 
#£(F P ) tells us that 

r p (D)=p ^ #E(W p )=p. 

Thus in standard terminology, the primes p > 7 in A(D) are exactly 
the anomalous primes for the elliptic curve E. 
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3. Arrows in the Index Divisibility Graph 

This section contains our main results. In Theorem 15 we classify 
the arrows (n — > nd) G Arr(D) for a large class of EDS, as described in 
the following definition. 

Definition. Let D be a minimal EDS associated to the elliptic curve 
and point (E,P). We say that D is 2-irregular if the following five 
irregularity conditions are true: 

(Ii) E has good reduction at 2, (I 2 ) #£(F 2 ) = 4, (I 3 ) r 2 = 4, 
(I 4 ) D 2 is odd, (I 5 ) ord 2 (L> 4 ) = 1, 

If any of the conditions (Ii) (I 5 ) is false, then we say that D is 2-regular. 
If in addition we have 

P G E ns (¥ p ) for all primes p | Disc(-E), 

then we simply say that D is regular. 

Remark 13. Our main result, Theorem 15, gives a good description 
of the index divisibility graph S(D) for regular EDS. Our decision to 
restrict attention to regular EDS represents a compromise between our 
desires for generality and conciseness, as well as the need to keep our 
exposition to a reasonable length. We remark that much of our anal- 
ysis goes through for non-regular EDS, in the sense that Theorem 15 
is still true for many (but generally not all) values of n, and that a 
long case-by-case analysis would give a lengthy statement that applies 
to most (maybe even all) values of n. In any case, we note that ev- 
ery D = (-D„) n >i contains a regular subsequence D' = (D n k) n >i, and 
then Theorem 15 applies to this subsequence. 

We start with a description of the index divisibility set of an EDS 
that will be a key tool for our classification. Its proof uses only the 
formal group properties of an EDS (Lemma 5). 

Proposition 14. Let D be a minimal regular EDS associated to the 
elliptic curve and point (E,P). Then the following are equivalent: 

(a) n | D n , i.e, n G 5(D). 

(b) There is some exponent e > 1 such that n \ D e n . 

(c) Every prime dividing n also divides D n . 

(d) For all primes p, we have p \ n =^> r p \n. 

Proof. Statements (b) and (c) are obviously equivalent, and (c) and (d) 
are equivalent by Lemma 4. It is also clear that (a) implies (b). It 
remains to show that (b) implies (a), i.e., that 

n\D e n — > n | D n . 
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It suffices to prove that for all primes p we have 

p | gcd(n, D n ) =>• ordp(At) > ord p (n). (7) 

So we let p be a prime dividing both n and D n and we write n = p u k 
with p \ k and v > 1. If z/ = 1, then (7) is obviously true (note p \ D n ), 
so we may assume that v > 2. 

We consider first the case that p | D pk . Applying Lemma 5(a) to 
D n = D pV -i. pkl we obtain 

ord p (A0 = ord p (D rk ) > oid^p^Dpk) >v = ord p (n). 

This shows that (7) is true in this case. 

We next suppose that p \ D pk , and we will show that either (7) is 
true or else D is 2-irregular. The assumption that p \ D pk is equivalent 
to r p \ ph. But we are assuming that p \ D p v k , so we have r p \ p v k. 
It follows that p 2 | r p , which is a very strong condition. In particular, 
since the regularity assumption implies that that P G E ns (¥ p ), and 
since r p is the order of P in E(¥ p ), we find that 

p 2 | #£ ns (F p ). 

Hence E has nonsingular reduction modulo p, and using the Hasse- 
Weil estimate, we further deduce that p = 2 and r 2 = #£ , (F 2 ) = 4. 
This gives conditions (Ii), (I2), and (I3) in the definition of 2-irregula- 
rity. Further, D 2 must be odd, since otherwise r 2 would divide 2, so 
we get condition (I 4 ). 

Since 2 | D 4 , so 2 | D ik , we can apply Lemma 5(a) to D n = D 2 v~z. Ak 
to obtain 

ord 2 (D n ) = ord 2 (ZV fc ) > ord 2 (2"- 2 J D 4 fc) = ord 2 (D 4k ) + v - 2. 

If ord 2 (_D 4 ) > 2, then this implies that ord 2 (D n ) > ord 2 (n), so (7) is 
true and we are done. Otherwise ord 2 (-D 4 ) = 1 and we have verified 
condition (I 5 ) for D to be 2-irregular. This is a contradiction, since we 
have assumed that D is regular, which completes the proof of Proposi- 
tion 14. □ 

We are now ready to state and prove our main theorem. 

Theorem 15. Let D be a minimal regular EDS associated to the elliptic 
curve and point (E, P). 

(a) Let n > 1. Then 

n E S(D) and ( P ,[!? n ° T f ) =^ (n np) E Arr(D). 
v ' \additive reduction J v ' v ' 
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(b) Let n > 1 and d > 1. Then 

neS(D) and d G A^(D) =► (^(/^L T t d) " 

Furthermore, 
ti G 5(D) and d G «4 ge n(D) and gcd(d,n) = l 

=>• (n ->■ nd) G Arr(D). 

(c) Le£ n > 1 and /e£ p be a prime such that 

n e 5(D), pfAi, arec ^ (n — >■ np) G Arr(D). 

(1) If E has good reduction at p and #E(F P ) ^ 2p, then 

P eA(D). 

(If P > 7, then we always have #E(F P ) ^ 2p.) 

(2) If E has bad reduction at p, then 

E has additive reduction at p. 

(d) Let n > 1 and d > 1 with d composite. Define 

t = (number of primes p \ d such that r p is composite), 
and 

Po = (smallest prime divisor of nd). 
Suppose that 

(n ->■ nd) G Arr(D). 

Then one of the following statements is true: 

(i) t = and d G „4 gen (D). 

(ii) t > 1 and 



n(.4)^n*f^. « 

Proof, (a) Suppose first that n G 5(D) and p \ D n . Write n = p l k with 
p\k. Then 

ordp(_D np ) > ord p (-D n ) + 1 from Lemma 5(a), 

> ordp(n) + 1 since n G 5(D), i.e., n \ D n , 

= i + l. 

Further, k \ n \ D n \ D np . Hence p l+l k \ D np , i.e., np \ D np , so 
np G 5(D). And since there are no proper divisors between n and np, 
it follows that the directed graph 5(D) contains the arrow n — > np. 
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Next we consider the case that n G 5(D) and p \ D n and E has 
additive reduction. Additive reduction implies that #E ns (F p ) = p, so 
that r p | p and p \ D p | D np . Meanwhile, n | D n \ D np by assumption. 
Since p \ D n , it must be that p \ n. Hence np \ D np and np G 5(D), 
from which we conclude that (n — > np) G Arr(D). 

(b) Let d — p\ • • • pt G -A gC n(D) be a generalized aliquot number for D. 
We will show that nd G 5(D). First, let p = pi be one of the primes 
dividing d. The rank of apparition satisfies r p . = Pi + ±, where for nota- 
tional convenience we let p^+i = p 1 . Hence Pi \ D Pi+l \ D nd . Next let p 
be a prime dividing n. Then p\n \ D n \ D nd . We have shown that any 
prime p dividing nd satisfies p \ D nd . By Proposition 14, we conclude 
that nd G 5(D). 

Now we determine the possible values of gcd(<i, n). If pi \ n, then 
since n G 5(D), Proposition 14 implies that pi + \ = r Pi must divide n. 
Hence, by the construction of d, we see that either n is divisible by 
none of the prime dividing d, or it is divisible by all of them. Therefore 
gcd(<i, n) — 1 or d. 

Suppose now that gcd(d, n) = 1 and that e | d is a divisor of d such 
that ne G 5(D). Then by the reasoning of the last paragraph, with n 
replaced by ne, we find that gcd(<i, ne) = 1 or d. But gcd(<i, n) = 1 
and e | d, so we conclude that e = 1 or e = d. Hence (n — > nd) G Arr(D) 
by definition. 

This completes the proof of (b). We also note that some condition 
such as gcd(n, d) — 1 is necessary. For example, suppose that (p, q) is 
an amicable pair and that p divides D n . Then there is no arrow from 
n to npq, because there are "shorter" arrows n — > np — > npq. 

(c) We are given that n G 5(D), np G 5(D), and p \ D n . Since 
np G 5(D), Proposition 14 implies p \ D np . We observe that 

p\D n <^=^ [n]P ^ O (mod p), 
p | D np <^=^> [np]P = O (mod p). 

Hence under our assumptions, in particular the regularity assumption, 
we see that the point [n]P has order exactly pinPe E ns (¥ p ). Hence 
p\E QS {¥ p ). 

(c-1) Suppose first that E has good reduction at p, so i? ns (Fp) = E(¥ p ). 
We want to show that p G ^4(D). The assumption that #E(F p ) ^ 2p, 
combined with Hasse's estimate \p + 1 — #E(F P )\ < 2^/p, implies that 

p\#E(W p ) ^ # E (W p )=p. (9) 

Since r p \ #E(F P ), we see that r p = 1 or r p = p. But r p = 1 implies 
that p | D 1 , contradicting p \ D n . Therefore r p = p, which implies that 
p G A(D), i.e., (p) is an aliquot cycle of length one. 
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(c-2) Next suppose that E has bad reduction at p. It follows from p \ 
Ens(Fp) that E has additive reduction at p. (If it had multiplicative re- 
duction, then E ns (¥ p ) would contain p±l points, depending on whether 
the reduction is split or nonsplit.) 

(d) We first show that gcd(-D n , d) = 1, which in particular implies 
that gcd(n, d) = 1, since n \ D n . To see this, suppose to the contrary 
that gcd(D n , d) > 1, and let p be a prime dividing gcd(-D n ,<i). Since 
p | D n , we know from (a) that (n — > np) G Arr(D). But since p \ d, 
we have divisibilities n \ np \ nd, so the fact that n — >■ np and n — )• nd 
are arrows implies that either n = np or np = nd. Neither of these is 
possible, since p > 2, and d is composite by assumption. This completes 
the proof that gcd(D n , d) — 1. 

In order to analyze the arrow (n — > nd), we associate to the integer d 
a directed graph Qd as in the following lemma. The graph Qd classifies 
the primes dividing each rank of apparition r p . 

Lemma 16. Let D be a minimal regular EDS, let n > 1 and d > 1 

with d composite, and assume that (n — > nd) G Arr(D). We construct 
a directed graph Q d with vertices and arrows defined as follows: 

Ver(Q d ) = {primes p such that p \ d}, 
Arr(^) = {p — > q : q \ d and q \ r p }. 

(N.B., the graph Qd is entirely distinct from the graph on S(D).) 

(a) Every vertex of Qd has an in-arrow. 

(b) Every vertex ofQd has an out-arrow. 

(c) The graph Q d is connected. 

Proof. Let q \ d be a prime divisor of d, i.e., q is a vertex of Qd- 

(a) We need to show that q has an in-arrow. Let d! = d/q. Since 
(n nd) G Arr(D), we know that nd' ^ <S(D). By Proposition 14, 
this implies the existence of a prime p | nd' satisfying r p \ nd' . Since 
r p | nd by Proposition 14, this implies that q | r p , which shows p — >■ q 
as required. 

(b) We need to show that q has an out-arrow. Since q | d, we have 
q | D d or r q \ d. This shows that some prime p \ d satisfies p \ r q , and 
thus q — > p as required. 

(c) Define d! to be the part of d supported on primes appearing as 
vertices in a connected component Q' of Qd- Then for each prime p \ d' , 
all primes dividing r p appear in Q' by connectedness. Since r p \ nd by 
the assumption that nd G 5(D), this implies r p \ nd' . This shows that 
nd' G 5(D), which contradicts {n — > nd) G Arr(D) unless d! — 1 or 
d! = d. So Qd is connected. □ 
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Suppose first that r p is prime for every p \ d. We need to prove 
that d G A gcn (D). By definition, for every arrow (p — >■ ?) € Arr(C? d ) we 
have q \ r p , so the assumption that r p is prime implies that r p = q. In 
particular, every vertex in the finite directed graph Qd has at most one 
outgoing arrow. But Lemma 16(b) tells us that every vertex in Qd has 
at least one outgoing arrow, and Lemma 16(c) says that the graph is 
connected. It follows that Qd consists of a single loop, 

Pi -> P2 -> P3 -> > Pt 




This loop satisfies r Pi = Pi+i, so by definition (pi, . . . ,p t ) is a general- 
ized aliquot cycle for D, and hence pip 2 • • -pt G ^4 ge n(D). Since we also 
know that gcd(<i, n) = 1, it follows from part (c) of the theorem that 

(n — > npxpi ■ ■ - pt) G Arr(D). 

But npip 2 ---pt | rirf, so that fact that (n — > nd) G Arr(D) implies 
that d = pip 2 • • - Pt- Hence d G ^4 gcn (D). This completes the proof of 
part (i). 

In order to analyze the case that one or more of the r p are composite, 
for each vertex q G ^ we let 

InDeg(g) = #{p | d : (p -> q) G Arr(^ d )} 

denote the in-degree of q, i.e., the number of arrows pointing in to q; 
and similarly OutDeg(g) will denote the out-degree of q. Lemma 16 
tells us that InDeg(g) > 1 for all q G Gd- For each p G Qd we know 
that r p is divisible by the primes at the tips of the outgoing arrows 
from p, so we can factor r p as 

r p = ( ] [ g^Mp for some M p > 1. 

Further, from Proposition 14, the fact that nd G 5(D) and p | d implies 
that r p I nd, so every prime divisor of M p is also a prime divisor of nd. 

We now multiply over all p G C^, i-e., over all p \ d, and rearrange 
the terms to deduce that 

n^=n( n M=(u^ m M m ')- 

p\d p\d q I d such that q\d p\d 

(p^g)eArr(CJ d ) 

Since InDeg(g) > 1 for every q \ d, we can rewrite this as 

Uj = {U^ Dcg{9) ~ 1 ){Il M ^ ( 10 ) 

p\d q\d p\d 
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where the right-hand side is a positive integer. Using the Hasse-Weil 
bound (^/p+1) 2 > #E ns (Wp) and the fact that r p | #£ ns (F p ), we obtain 
the useful inequalities 

n l + 1 v > n > (nr*'-) (n*4 m 

p\d X VF/ p\d F q\d p\d 

We now use (11) to derive a bound that depends on the number 
of composite r p . (See also Remark 19.) Let p be the smallest prime 
divisor of nd. Then 

JJ ? InDeg( ? )-l > -Q p InDcg(g)-l = ^E,| d (InDcg( g )-l) 

<zM i\d _ #KTT(g d )-#v C v(g d ) 

— Po 

= p E p|d (OutDe g(p) -l)_ (12) 

Now consider a prime p & Qd such that r p is composite. If r p is divisible 
by two or more primes that also divide d, then OutDeg(p) > 2, so we 
get a factor of p in (12). On the other hand, if there is some q | r v 
with q \ d, then q | M p , so we get a factor of g in (11). Further, we 
must have q \ n, since as noted earlier, r p | nd. Thus q > Po- This 
proves that every composite r p with p | d contributes a factor to (11) 
that is greater than or equal to p . Hence the lower bound in (11) is at 
least pi, where t is the number of p \ d such that r p is composite. □ 

The following corollary may be compared with Smyth's result [19, 
Corollary 2] for Lucas sequences. 

Corollary 17. Let D be a minimal regular EDS, let n e <S(D), and 

let m be an integer of the form 

m = pip 2 •■•Ps- did 2 ■••d t , 
where the primes Pi and integers di satisfy 

Pi | D n and d { e w4 gen (D). 

Then nm e 5(D). 

Proof. This is immediate from Theorem 15(a,b) and induction on the 
number of factors of m. □ 

4. Remarks on Arrow Construction 

In this section we make a number of remarks concerning the existence 
of index divisibility arrows as described in Theorem 15, and we give 
examples of non-standard arrows as per Theorem 15(d-ii). We assume 
throughout that our EDS is minimal and regular. 
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Remark 18. Given an element n G 5(D), Theorem 15 gives two "stan- 
dard" ways to create arrows (n — > nd) G Arr(D). First, Theorem 15(a) 
gives an arrow (n — > np) for each prime p \ D n . Second, Theorem 15(b) 
gives an arrow for each aliquot number d G ^4 gen (D) that is prime to n. 
Conversely, Theorem 15(d) implies that any "non-standard" arrow sat- 
isfies 



n 

P \d 



(13) 



In particular, writing v(d) for the number of distinct prime divisors 
of d and p m i n {d) for the smallest prime dividing d, we have 

| log 2 log 2 



VPmin(rf) +0(1). 



log(l+p min (d)-V2) 2 

Thus if the smallest prime divisor of d is large, then u(d) will be large, 
and d will be enormous. The following brief table uses (13) to give the 
smallest values of v{d) and d for various values of p m m(d). 



Pmin(d) > 


10 


10 2 


10 3 


10 4 


10 5 


u{d) > 


2 


4 


12 


36 


100 


d> 


143 


1.21 • 10 8 


1.56 • 10 ae 


1.80 • 10 i44 


1.85 • 10 5UU 



And if Theorem 15(d) gives a lower bound for (13) that is larger 
than v^2, then the lower bounds for v(d) and d in terms of p m - m (d) 
will be even larger. 

Remark 19. The formulas (10) and (11) derived during the course of 
proving Theorem 15(d) impose stringent conditions on the allowable 
values of d. We used these formulas to derive a general lower bound, 
but when analyzing a specific EDS, it is probably best to use them 
directly. We also note, although we will not prove, that (10) is true 
even if d is divisible by primes for which P has singular reduction. 
Similarly, the following version of (11) is true in general: 

n^^n^-'Xn^). 

p\d " q\d p\d 

where € is the Neron model of E. Note that if E has bad reduction, 
then #(£(F P ) = c p #E QS (F p ), where c p is the number of components in 
the special fiber above p. In particular, c p < 4 unless the reduction is 
split multiplicative, in which case c p = ord p (Disc(i?)). 

Example 20. Continuing with the EDS associated to the elliptic curve 
and point from Example 1, we have 

Disc(E) = 37, #£(F 2 ) = 5, #£(F 3 ) = 7, #£(F 5 ) = 8. 
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In particular, E has multiplicative reduction at 37 and good reduction 
elsewhere, the point P is in i? ns (F 37 ), and #E(W P ) ^ 2p for all primes p. 
Further, since D 5 = 2 and D w = 4, we see that Lemma 5(b) is true 
even for p = 2 and all values of n and A;, so we can treat 2 as we do all 
other primes. 

We claim that for all primes p, 

p\D n or p e A(D) ^=> (n np) G Arr(D). 

The implication =>- follows directly from Theorem 15(a,c). Conversely, 
if (n — > np) G Arr(D), then either p \ D n , or else Theorem 15(c) tells 
us that p E A(D). We thus have a precise description of the arrows of 
prime weight. 

Theorem 15(d) says that arrows {n — > nd) of composite weight with 
d ^4 gcn (D) have d values that are either divisible by small primes 
or are huge. Further, examining the proof of Theorem 15(d) shows 
that the prime divisors of such d must satisfy some fairly stringent 
conditions. We suspect that for this example there are no such arrows, 
i.e., 

or d G «4 g en(D) <=> (n nd) G Arr(D). 

Example 21. The following example shows that "non-standard" arrows 
exist (cf. Remark 18). Let D be the EDS associated to 

E : y 2 + 2xy + y = x 3 + x 2 + 7x + 4 and P = (4,7). 

The curve E is nonsingular at 2, 3, and 5, and 

#£(F 2 ) = 3, #£(F 3 ) = 5, #£(F 5 ) = 6. 

Further, the point P has exact order 6 in E(¥ 5 ). Thus r 2 = 3, r 3 = 5, 
and r 5 = 6, so 

2,3,5,6,10,15 ^ 5(D) and 1,30 G 5(D). (14) 

Alternatively, we can verify (14) directly by explicitly computing the 
relevant terms of D, 

Di mod 1 = 0, D 2 mod 2 = 1, D 3 mod 3 = 2, D 5 mod 5 = 4, 
D 6 mod 6 = 4, D w mod 10 = 3, D 15 mod 15 = 3, -D 30 mod 30 = 0. 

It follows from the definition of the directed graph 5(D) that (1 — > 
30) G Arr(D). However, since r 5 = 6 is not prime, we have 30 ^ 
A gen (D). Thus the arrow (1 — > 30) is not predicted by Theorem 15 (d- 
i). This does not contradict the theorem, of course, since 

#£(F 2 ) #£(F 3 ) #£(F 5 ) _ 3 5 6 
2 ' 3 ' 5 2 ' 3 ' 5 



d is prime \ 
and d | D n l 
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so condition (8) is satisfied and we are in the situation of Theorem 15 (d- 
ii). 

Remark 22. Generalizing Example 21, we sketch how to construct 
EDS having non-standard arrows with arbitrarily large values of d. 
The proof of Theorem 15(d) suggests the method. We start with 
primes pi, . . . ,p N and integers ni, . . . , and ki, . . . , k N satisfying 

\pi + l - kiUi\ < 2y/p~. 

Our goal is to find an elliptic curve E/Q and point P G E(Q) such 
that #E(F Pi ) = kiiii and r Pi = rii for all 1 < i < N . 

A theorem of Deuring [4] says that there exists an elliptic curve 
Ei/¥ Pi satisfying 

# J E? i (F p .) = km, 

and a result of Ruck [15, Theorem 3] says that we can choose Pj so 
that the group structure of Pj(F p .) ensures the existence of a point p G 
Ei(¥ p ) of order rii. Making a change of coordinates, we may assume 
that Pi = (0,0). 

Next we apply the Chinese remainder theorem to the coefficients of 
the Weierstrass equations of E\, . . . ,E n . This gives an elliptic curve 
E/Q with (0,0) G E(Q) that satisfies 

E mod Pi = E t , 1 < i < N. 

If the Weierstrass equation for E is not globally minimal, then we can 
change coordinates to make it minimal without affecting the reduction 
at pi, . . . ,pjv, since they are primes of good reduction. For simplicity, 
we will assume that some n« is divisible by a prime greater than 7, 
since then Mazur's Theorem [16, VIII. 7.5] ensures that (0,0) is not a 
torsion point. We may thus associate to E and P an elliptic divisibility 
sequence D = (D n ) n >i satisfying 

r p . = rii for all 1 < % < N. 

Finally, we observe that arbitrarily large non-standard arrows can 
be constructed in this way. We begin with any prime p±, we let 

p u p 2 ,...,p N 

be a list of consecutive primes, and we set 

d = pjp 2 p3 ---Pn- 

We then find a curve and point whose associated EDS satisfies 

r Pi =p i+1 , l<i<N -1, and r PN = p\. 

If the list of primes is taken to be long enough, then the final con- 
dition r PN = p\ is allowed by Hasse's bound, and we can proceed 
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as in the description above to find a sequence D = (D n ) n >i with 
(1 -> d) G Arr(D). 

Example 23. We use the method described in Remark 22 to construct 
a non-standard arrow (1 — >■ d) for the moderately large integer 

d = 5 2 - 7- 11 -17 = 32725. 

We want to construct an elliptic curve E/Q and point P G E(Q) 
satisfying 

r 5 = 7, r 7 = 11, m = 17, and r i7 = 25. (15) 

Then the associated sequence D = (D n ) n >i will have (1 — > d) G Arr(D), 
according to Proposition 14. 

To do this, we first found elliptic curves E 5 /¥ 5 , E 7 /¥ 7 , En/Wu and 
En/Fn satisfying 

#£ 5 (F 5 ) = 7, #£ 7 (F 7 ) = 11, #Sn(F n ) = 17, #£ 17 (F 17 ) = 25. 

This is possible because the Hasse bound is satisfied in each instance. 
We then used the Chinese remainder theorem to find an elliptic curve 
E with minimal Weierstrass equation 

y 2 + y = x 3 + x 2 - 1291874622406186a; + 17872226251073822113702, 
and point 

P= (20751503,1073344). 

(We've moved P away from (0, 0) to make the numbers a bit smaller.) 
The associated sequence D = (D n ) n >i begins 

1, 2146689, 286883381041833542301, 

60768120452650698495048133538894517, . . . 

By construction, (1 ->■ 32725) G Arr(D). Of course, the 32725 th term is 
too large to print, but the claim can be verified by computation modulo 
32725. 

For this example, we can verify equation (8) in Theorem 15(d), which 
states 

In our case, p — 5, t — 1, and the left-hand side exceeds 10. 

Remark 24. In the definition of EDS, the elliptic curve may be re- 
placed with a singular cubic curve as long as P is a non-singular point, 
since E ns (Q) is a group. More precisely, E ns (Q) is either the additive 
group Q + , the multiplicative group Q*, or a subgroup of a quadratic 



Index divisibility in elliptic divisibility sequences 



21 



twist of the multiplicative group; see [16, III. 2. 5, Exercise 3.5]. Thus 
EDS on singular elliptic curves are closely related to Lucas sequences. 
For example, consider the nodal singular cubic curve and point 

C : y 2 + 3xy + 3y = x 3 + 2x 2 + x, and P = (0, 0). 

The associated EDS, 

D : 1, 3, 8, 21, 55, 144, 377, 987, 2584, 6765, . . . 

consists of the even-indexed Fibonacci numbers. This is exactly the 
Lucas sequence generated by 

L n +2 — 3L n+ i — L n , Lq = 0, L\ — 1. 

The index divisibility set of D is 

«S(D) = {1, 5, 6, 12, 18, 24, 25, 30, 36, 48, 54, 55, 60, 72, 84, . . .}. 

In the notation of Smyth's Theorem 2, we have 

a = 3, 6 = 1, A = 5, and B 3A = {1 ->■ 6}. 

In the language of our paper, 5, 6 G ^4 gcn (D), since 

r-i = 3, r 3 = 2, and r 5 = 5. 

Thus (2, 3) and (5) are generalized aliquot cycles. Notice that the curve 
C reduces modulo p to a curve having p, p — 1 or p + 1 non-singular 
points according as p ramifies, splits, or is inert in Q(y / 5)- 

In general, our Theorem 15 and Smyth's Theorem 2 can probably be 
combined into a general theorem on (possibly singular) cubic curves. 
Notice that Smyth's set B a ,b m &y include non-standard arrows in the 
case of the multiplicative group, although the analysis is simpler be- 
cause #C ns (F p ) e {p,p+ l,p — 1}- The primes p dividing A = a 2 — 46 
are the primes for which the group underlying the Lucas sequence re- 
duces to the additive group F+. They are thus analogous to the primes 
of additive reduction whose arrows (n —> np) are described in The- 
orem 15(a,c). We also note that in the multiplicative group case we 
never have r 2 = 4, so we are always in the 2-regular setting. 

5. Elliptic aliquot cycles 

Let D be an EDS with associated elliptic curve and point (E,P), 
and let (p, q) G A(D) be an amicable pair for D. Then the point P has 
order q modulo p, and P has order p modulo q. This implies that 

q | #E(F p ) and p | #E(W q ). 
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Conversely, if we are given D and (E,P), and if p and q are distinct 
primes of good reduction satisfying 

#E(W p ) = q and #E(W q ) = p, (17) 

then (p, q) is automatically an amicable pair for D. 

We note that the conditions (17) do not refer to the point P. This 
leads to the following definitions. 

Definition. Let E/Q be an elliptic curve. An aliquot cycle of length £ 
for E/Q is a sequence (pi,P2, • • • ,Pe) of distinct primes such that E has 
good reduction at every pi and 

#E(W Pl )=p 2 , #E(W P2 )=p 3 , #E{W Pl _ 1 )=p t , #E{W Pt )=p 1 . 

An amicable pair for E/Q is an aliquot cycle of length 2. 

Remark 25. The distribution of amicable pairs and aliquot cycles on 
elliptic curves is studied in [17]. In particular, it turns out that elliptic 
curves with complex multiplication behave quite differently from curves 
without CM. For the convenience of the reader, we briefly summarize 
some of the material in [17]. 

• If E(Q) contains a non-trivial torsion point, then E has (essen- 
tially) no aliquot cycles. This is clear since E(Q) tOTS E(F p ) for 
all primes p \ 2 Disc^/^; cf. [17, Remark 5]. 

• For any £, there exists an elliptic curve E/Q that has an aliquot 
cycle of length £. More generally, for any £i, . . . ,£ s there exists an 
elliptic curve having disjoint aliquot cycles of length £i, . . . ,£ s [17, 
Theorem 13]. 

• Let E/Q be an elliptic curve with complex multiplication and 
j(E) 7^ 0. Then E has no aliquot cycles of length £ > 3 composed 
of primes p > 5 [17, Corollary 16]. 

• Let E/Q be an elliptic curve with j(E) = 0. Then E has no aliquot 
cycles of length 3 composed of primes p > 11 [17, Proposition 48]. 

• Conjecture: Assume that there are infinitely many primes p such 
that #E(F p ) is prime. If E does not have CM, then 

#{aliquot cycles (pi, ...,pe) with p, < X} »< 

(log A J* 

If E has CM, then there is a constant Ce > such that 

X 



^{amicable pairs (p, q) with p,q < X} ~ Ce- 



'(log*) 2 

The next proposition shows that aliquot cycles for an elliptic divis- 
ibility sequence are closely related to aliquot cycles on the associated 
elliptic curve. 
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Proposition 26. Let D be a minimal EDS, and let (E, P) be the as- 
sociated elliptic curve E/Q and point P e E(Q). 

(a) Let (pi, . . . ,pe) be an aliquot cycle for E/Q such that pi \ D 1 for 
all i. Then (pi, . . . ,pi) is an aliquot cycle for D. 

(b) Let (pi, . . . ,pe) be an aliquot cycle for D. Then 

T]i^M<2 =► (fy--y p <\ i8 , an tP ). (is) 

J- J- p i yahquot cycle for h J v ' 

In particular, 

min v . > \ =^ ( -'ft)" " ^ fig) 

i™/' (2V2< _ 1)2 ^ ^a%uo£ cyc/e for EJ ' { ' 

cf. Theorem 15(d). 

Proof, (a) If (pi,p2, ■ ■ ■ iPe) is an aliquot cycle for -E/Q, then for all % 
we know that #E(F Pi ) = pi + \ is prime. Since pj + i \ D±, the order of 
the point P in E(F Pi ) must equal Pi+i. Therefore r Pi (D) = Pi+i, so the 
cycle is aliquot for D. 

(b) The proof is similar to the proof of Theorem 15(d). We are given 
that r Pi (D) = p i+ i for all i, or equivalently, the point P has order p i+1 
in the group E(¥ p .). Thus for every 1 < % < I we have 

#E(F Pi ) = p i+1 M Pi for some M Pi > 1. 

Multiplying for 1 < i < t and dividing by p\ ■ ■ • pn yields 

i=i p * »=i 

Thus the assumption that f|i #E(¥ Pi )/pi < 2 implies that Mj = 1 for 
every i, so (pi, ■ ■ ■ ,pe) is an aliquot cycle for E 1 . This proves (18). 

To prove (19), we use the Hasse-Weil bound #E(W P ) < (y/p + l) 2 to 
obtain 



Now a little bit of algebra, combined with (18) yields (19). □ 

6. Miscellaneous Remarks 

We conclude with two brief remarks. 

Remark 27. Recall that a sequence A = (A n ) n >i is called a divisibility 
sequence if 

m n — V A m A n . 



< 1 + 



24 



JOSEPH H. SILVERMAN AND KATHERINE E. STANCE 



Examples of divisibility sequences include Lucas sequences of the first 
kind, the odd terms of Lucas sequences of the second kind, and elliptic 
divisibility sequences. We observe that if A is a divisibility sequence, 
then 

n e 5(A) and d | D n and gcd(n, d) — 1 =>- nd E S(A). 
In particular, there is a sequence of arrows in Arr(A) satisfying 

n — > ■ ■ ■ — > nd. 

This is one way in which the index divisibility graph of divisibility 
sequences exhibits a structure not found for arbitrary sequences. It 
might be interesting to see if there are any other general statements that 
one can make about the index divisibility graph of general divisibility 
sequences. 

Remark 28. A classical alternative definition of an elliptic divisibility 
sequence is a sequence of integers W = (W n ) n >i defined by four initial 
terms {W\, W 2: W 3 , W 4 ) and satisfying the recursion 

W n+m W n - m W* = W n+r W n - r W^ - W m+r W m _ r Wl for all n > m > r. 

One can show that if the sequence is normalized by W\ = 1 and 
W 2 | W4, then every term is an integer. Ward [25, 26] was the first 
to study the arithmetic properties of these sequences. Subject to 
some non-degeneracy conditions, he showed that there is an elliptic 
curve E/Q given by a Weierstrass equation and a point P E E(Q) 
such that 

w n = MP), 

where ip n is the n'th division polynomial for E [16, Exercise 3.7]. (See 
[25] or [18, Appendix A] for explicit formulas for E and P in terms of 
the initial terms of the EDS.) In particular, if D = (D n ) n >i is the EDS 
associated to (E, P), then D n \ W n for all n > 1. Thus 

n I D n n I W n , (20) 

so index divisibility for D is a stronger condition than it is for W. 
Further, one can show that 

ord p (L> n ) = ord p {W n ) 

for all primes p at which the Weierstrass equation has good reduction, 
so the implication (20) can be reversed if we ignore primes of bad 
reduction. This shows that the divisibility properties of D and W are 
closely related. We have chosen in this paper to concentrate on the 
former. 
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